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Introduction 


¢ Providers of externally-mounted scientific payloads at the 
International Space Station (ISS) are required to evaluate column 
number density (CND, o) associated with various gas releases and 
demonstrate that they fall below some maximum requirement 
— Must be considerate of other payloads 
— Since this includes unknown future additions, becomes a search for 
maximum CND along any path 
¢ Occasionally astrophysicists are interested in estimating the 
amount of gas released by some event or process by evaluating 
light attenuation of a distant star having known properties due to 
this release 
— Milky Way center, black hole, “Fermi Bubbles” 
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ESA/Planck Collaboration (microwave); NASA/DOE/Fermi LAT/Dobler et al./Su et al. (gamma rays) 
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¢ Develop analytical CND expressions for general paths that 
intercept various common point sources under high vacuum 
conditions 
— Effusion/low rate evaporation/outgassing (M = 0) 
— Venting via sonic orifice (M = 1) 


— Spherically-symmetric, radial expansion (M = 0) 
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“Venting (Directed) Source Behavior 





¢ External neutral gas phase sources on ISS result from a 
number of different physical mechanisms 
— Supersonic expansion through thruster nozzles 
— Pressure-driven acceleration to sonic conditions across an orifice 
— Surface evaporation, desorption (may or may not have bulk velocity) 
— Effusion--low-rate, high-Kn venting (M = 0) 
— Diffusion-limited outgassing (M = 0) 
¢ This study assumes that, for these applications, the point 
source may be described using free molecule flow model 
approximations 
— Density levels fall rapidly with distance from source location 


— Existence of self-scattering collisions may not substantially alter plume 
distribution from free molecule flow description 
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Directed Source—Steady Density 





¢ Can compute many different types of local quantities at receiver 
position x relative to source 


— Steady number density n from a directed axisymmetric source given by 


a 9 2 
n(x,1)= POS? ow me je + S40? a (vet | 
Ar 2 
Ue 


— Speed ratio s = Gu, aa scos@ 


e 


— A,: normalization factor, function of s 
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“Column Number Density (CND, o) 





¢ Integrated effect of molecules encountered across a prescribed 
path / 
— When unbounded, 
CO 
o =([ndl 
0 
¢ For ISS application, the requirement not to exceed 0,,; allows 
one to determine the physical envelope around the source 
where the limit is violated 


— With a singularity at the source origin, the model will always predict 
some critical envelope 


¢ Not consequential for low N 
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— Effusive CND Expressions 


¢ For low rate, high-Kn venting through an orifice with thermal 


effusion, no bulk motion, plume model density simplifies to 
ne 0) = N cos @ 

2 
r~ V87RT 


— Also describes density field due to outgassing or low rate volatile 
evaporation from a planar surface viewed from a distance 


¢ Column number density given by 


N &cos@ 


——— A dl 
SART 0 r 





Oo= 
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1-D Centerline Path 


l=r—-Xxo 






Source surface plane 


N 


XO / S7zRT 


* For effusion, the centerline result is simply o,) , ( Xo ) = 


¢ Since density is maximized along centerline, tempting to consider this path 
produces the highest CND. However, this is not so! 
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““ 2-D Path, Surface Plane Intersection 





/sinn =rcos@ x 
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¢ Solution for effusion becomes 


N t 
fe siny a an 77 sin7 


V82RT 1% (1 — cos 1) es l1—cos7 


¢ In the limit where 7 > ©, 7 > 0 


— Vanishingly small distortion of triangle to describe a path parallel to 
source plane at height x, find 


o(n >0)> 2O cle 


— Special case may be confirmed by evaluating o along horizontal path at 
height x, directly 


— This case provides the maximum CND for effusion 
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rcos@=Isin@+ 7p cos & 





Initial location above a 
source surface plane we 
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¢ Effusive gas solution: 


N  sin@—cos 
o = —— 


V82RT 1% (1 — cos 77) 


¢ Solution still maximized for distant points along paths parallel 
to source plane separated by x, 


— Collapses to previous solution 


O. 


max,e >2 Ocle 
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¢ When bulk fluid motion is involved (s > 0), plume model 
behavior becomes too complex to handle directly (w = s cos6) 


n(x,t)= PE ee pe 1 Sew? a +e | 


Ayr? 


¢ Decided to approximate the model behavior, replacing angular 
distribution by cos’ 6 


cos? 0 


2 
r 





Ny (r, 0) ~K 
— Good approximation for many species, different types 
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Sonic Angular Distribution Comparison 


Svo gals 


——= model, gamma = 5/3 
——= model, gamma = 7/5 


——— model, gamma = 4/3 


== model, gamma = 9/7 


—m™ + cosr3 


04 - = — cos4 
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0.2 y=7/5, 4/3, & 9/7 





NORMALIZED RADIAL DIMENSION 


! 


0 0.2 0.4 0.6 0.8 1 
NORMALIZED AXIAL DIMENSION 
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: K 
¢ 1-D centerline case: o,),, =— 
XO 





e 2-D, M centerline & source surface plane: 


O, = = 20 + cos 1)+ z sin7 sn 
3% Sinn Z 





_ Ocl,s 
3cos7 





1 
20 + cos 1)+ 5 sin77 sin 2n| 


Source surface plane, 
edge on 





4 
— Maximum effect: o(n > 0) =F 3 Ochs 
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“= 3-D Path, Sonic Approximation 


¢ Generally, 


O¢j.s tan in? | —cos? | 
c= fast Ce —3sin? @cos Ay + 3sin@cos” ayes" (2-0) 


3(1 — cos n) (1 + cos n) 
¢ Maximum effect when 


Omax,s > 3 Ocl,s 


@ in plane // x 


s) 7in plane containing r, ry, & 1 





Initial location above so 
source surface plane awe 
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— Higher VM CND Observations 


~ cos” O 
- Assume adequate fit for our purposes using n (7,0) * K — 
K r 
¢ For axial, centerline case, find o,, =— 
x0 


¢ From previous results, might think limiting transverse case becomes 
" m4+1 
O xverse— Oc] 

Mm 


— Always larger than axial 


Actually oo r 
m+ 
O xverse — vib 5 mt). V10-) ——~ 





— Axial case is larger for m>5 
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Radial Point Source 


¢ Model spherically-symmetric expansion 
— No directional constraints 
— No bulk velocity (thermal expansion, s = 0) 
¢ Use solution due to Narasimha 
N 


- mr V 82RT 


ne 
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” Radial Point Source CND Expressions 


¢ Generally, 
TW Source 


N 
ry V8n0RT siny 


— Notice ry siny acts like x, in venting cases 





oOo= 





Yo 

¢ When w = z (path along source radial line) Initial location 
_ NN 

’ Thy 8ART 


¢ When yw = 7/2, path begins at right angles to source, 7) = Xp, and 


ype aes 
ees 5 Or 


¢ Maximum CND found for path that extends to infinity in both directions: 


Omax,r — 7 Oy (The m = 0 result!) 
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Concluding Remarks 


¢ Undertook a study to determine closed form analytical solutions for a number of 
frequently encountered CND configurations 


¢ For low-rate effusive venting and higher-rate sonic discharges, maximum CNDs 
should occur along paths parallel to the source plane that intersect the plume axis 
¢ Maximum CNDs for paths immersed in the presence of an unconstrained radial 
source do not lie along radial trajectories 
¢ For source angular distributions ~cos” @, it was shown for integer values of m>5, 
maximum CND values switched from transverse to axial paths 
— Likely associated with spacecraft thruster plumes 


¢ These analytical solutions and associated observations should greatly reduce the 
amount of effort needed to assess CNDs for a variety of space-related applications 
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Plume Model Formulation—Source 





¢ Find particular solution to collisionless Boltzmann equation 


for source Q;: 


Sy St, ge T - 
a Ox = 21 


where Q, represents a Lambertian source superimposed on a 
bulk velocity 


A 
QO, = “P_ 5(x)nle)y Alexp- B2(v — Ue ?) 
17% 


and the normalization factor is given by 


Dn 
Aj=e” °° bn scos¢.(1+ erf(scos@, )) 
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“Plume Model Formulation—Definitions 
2 









source Axis along exit normal # 


‘ receiver normal 
location ecole: donne 


receiver 
X,v location 





Subscript e represents exit conditions from source 


Simplifies for axisymmetric conditions 


?. = 0 
—~ ¢=0 
¢ other definitions: s = fu, = —"e _. w=scos0 
2RT. 


e 
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